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Abstract 

Heavy-hadron chiral perturbation theory (HHxPT) is applied to the decays of the even-parity 
charmed strange mesons, D s o(2317) and L> s i(2460). Heavy-quark spin symmetry predicts the 
branching fractions for the three electromagnetic decays of these states to the ground states D s 
and D* in terms of a single parameter. The resulting predictions for two of the branching fractions 
are significantly higher than current upper limits from the CLEO experiment. Leading corrections 
to the branching ratios from chiral loop diagrams and spin-symmetry violating operators in the 
HHxPT Lagrangian can naturally account for this discrepancy. Finally the proposal that the 
D s o(2317) (.D s i(2460)) is a hadronic bound state of a D (D*) meson and a kaon is considered. 
Leading order predictions for electromagnetic branching ratios in this molecular scenario are in 
very poor agreement with existing data. 
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I. INTRODUCTION 



The discovery of the D s0 (2317) [l] and -D s i(2460) Q has revived interest in excited 
charmed mesons. The dominant decay modes of these states are D s0 (2317) — > D s tt° and 



D s i(2460) — > D*tt°, with widths less than 7 MeV [2j. There is experimental evidence 
indicating that D s0 (2317) and -D s i(2460) are J p = + and 1 + states, respectively 3,0]. Had 
the masses of the + and 1 + states been above the threshold for the 5*-wave decay into D 
mesons and kaons, as anticipated in quark model [5|, |6( as well as lattice calculations ODD 
they would have had widths of a few hundred MeV. In reality, the unexpectedly low masses 
make those decays kinematically impossible. The only available strong decay modes violate 
isospin, accounting for the narrow widths. 

The _D s o(2317) and _D s i(2460) can also decay electromagnetically. The possible decays are 
L> s i(2460) -> £>*7, £> sl (2460) -> L> s 7, and A, (2317) -> D* sl . The decay D s0 (2317) -> D s7 
is forbidden by parity conservation. In the heavy-quark limit, both the three electromagnetic 
decays and the two strong decays are related by heavy-quark spin symmetry Q. Belle has 
observed the decay -D s i(2460) — ► D s ^ from D sl (2460) produced in the decays of B mesons |3|] 
and from continuum e + e _ production "\. The ratio of the electromagnetic branching fraction 
to the isospin violating one pion decay reported by the experiment is 



Br(D sl (2460) -> D s7 ) I 0.38 ± 0.11 ± 0.04 [3] 
Br(D sl (2460) - Dfifl) ~ j .55 ± 0.13 ± 0.08 [4]' 

In each case the first error is statistical and the second systematic. The other electromagnetic 
decays have not been observed. CLEO quotes the following bounds on the branching fraction 
ratios 

B,(IM2460) ^ P;,) Br( C , (2317) ^ C-oQ 

Br(£> s i(2460) -> L> s *tt°) Br(£> fl0 (2317) -> D s 7r°) v 7 

(The BELLE collaboration quotes weaker lower bounds of 0.31 and 0.18, respectively, for 
these ratios 

In this paper the decays of the _D s o(2317) and .D sl (2460) are analyzed using heavy-hadron 
chiral perturbation theory (HHxPT) HLLyPT is an effective theory applicable to the 
low energy strong and electromagnetic interactions of particles containing a heavy quark. 
It incorporates the approximate heavy-quark and chiral symmetries of QCD. Corrections 
to leading order predictions can be computed in an expansion in Aq CD /mQ, M/A x , and 



p/A x , where mg is the heavy quark mass, M is a Goldstone boson mass, p is the typical 
momentum in the decay, and A x is the chiral symmetry breaking scale. 

In section II, the leading order HHxPT predictions for the branching ratios are derived: 

1^(2460)^:7) - 37±0 07 Br(^ (2317) ^7) _ Q3 
Br(D sl (2460) - D*tt°) " ^ Br( J D s0 (2317) - D s tt°) " U ^ ' W 

Seading order calculations of strong and electromagnetic decays were first done in Refs. 
. Il4j].) These predictions deviate significantly from the CLEO limits. At next-to- leading 



order (NLO) there are 0(l/m<g) suppressed heavy-quark spin-symmetry violating operators 
as well as one-loop chiral corrections to the electromagnetic decays. Once these effects 
are included, predictions for the ratios in Eq. (J2J) can be made consistent with the present 
experimental bounds with coupling constants in the Lagrangian of natural size. 

The splitting between the even- and odd-parity doublets should be approximately the 
same for both bottom strange and charmed strange mesons. Therefore it is likely that the 
B s even-parity states will be below threshold for decay into kaons and narrow like their charm 
counterparts. The calculations of this paper can also be applied to the electromagnetic and 
strong decays of even-parity B s mesons when these states are discovered. 

The leading order HH;\PT Lagrangian used in this paper is invariant under nonlinearly 
realized chiral SU(3)l x SU(3)r and no further assumptions are made about the mechanism 
of chiral symmetry breaking. Models that treat the _D s o(2317) and D s i(2460) as the + and 
1 + chir al p artners of the ground state charm strange mesons are proposed in Refs. [13I 
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ia liana 

j. In these models, referred to as parity-doubling models, the even-parity and 
odd-parity mesons are placed in a linear representation of chiral SU(3)l x SU(3)r. These 
fields couple in a chirally invariant manner to a field S that transforms in the (3, 3) of 
SU(3)l x SU(3)r. The S field develops a vacuum expectation value that breaks the chiral 
symmetry. The resulting nonlinear sigma model of Goldstone bosons coupled to heavy 
mesons has the same operators as the HH%PT Lagrangian used in this paper. The assumed 
mechanism of chiral symmetry breaking in parity doubling models predicts relationships 
between coupling constants in the HH%PT Lagrangian. For example, the parity doubling 
models predict that the hyperfine splittings of the even- and odd-parity doublets are equal. 
This is in agreement with experimental observations. Other relationships between coupling 
constants in HHxPT are predicted Q, Q] by the theory of algebraic realizations of chiral 
symmetry (2^), in which hadrons are placed in reducible representations of SU (3) l x SU {3)r. 
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QCD sum rules have also been used to calculate some of the HH^PT couplings 21]. When 
more data on the electromagnetic decays of even-parity D s and B s mesons becomes available, 
the formulae derived in this paper can be used to extract the relevant couplings and test 
these theories. 
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the _D s o(2317) and P s i(2460) has prompted reexamination of quark 
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include DK molecules 
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as well as speculation that these states are exotic. Possibilities 
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35 1 . Masses have been calculated in lattice QCD 
(HQET) sum rules HQ, and potential as well as other models b 



, and tetraquarks |2 

heavy-quark effective theor 
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H H 0, 



The results of some of these papers are contradictory. For example, the lattice calculation of 
Ref. j3| yields a + -0~ mass splitting about 120 MeV greater than experimentally observed, 
quotes errors of about 50 MeV, and argues this is evidence for an exotic interpretation of 
the state. On the other hand, the lattice calculation of Ref. (stJ obtains similar numerical 
results but concludes that uncertainties in the calculation are large en oug h to be consistent 



with a conventional cs P-wave state. Some quark model analyses |2J 



26, 



4jJ conclude that 



interpreting the states as conventional cs P-wave mesons naturally fits the observed data, 



others reach the opposite conclusion (22 , 



42|- 



There have been some attempts to determine the nature of the P s o(2317) and P s i(2460) 



from the observed pattern of decays 



cays 



Q 44,0, 



3] as well as their production in 6-hadron de- 
l46j |. Ref. |23j] argues that the total width and electromagnetic branching 



ratios can distinguish between cs P-wave states and DK molecules, a nd g ives predictions 



l44j argue that the 



for these branching ratios calculated in the quark model, .nets. j^J 
observed branching fractions for B — > D s0 (2317)D^ and B — > P sl (2460)P^*^ are smaller 
then expected for cs P-wave states, suggesting that these states are exotic. These analyses 
assume an unproven (but plausible) factorization conjecture for the decays as well as quark 
model estimates for the P s o(2317) and D*i f 2460) decay constants, and have recently been 
extended to [0 and semileptonic B s decays 0]. 

Section III addresses the question of whether a model independent analysis of the decays 
can provide insight into the nature of the P s o(2317) and P s i(2460). In HHxPT the fields 
describing the + and 1 + mesons are added to the Lagrangian by hand. The only assumption 
made about these states is that the light degrees of freedom in the hadron are in the 3 of 
SU(3) and have j p = \ + . (In this paper, J p refers to the angular momentum and parity of a 



4 



heavy meson, and j p to the angular momentum and parity of the light degrees of freedom.) 
Light degrees of freedom with these quantum numbers could be an s quark in an orbital P- 
wave or sqq quarks all in an S"-wave. Therefore, a conventional quark model cs P-wave state 
and an unconventional csqq tetraquark state will be represented by fields having the same 
transformation properties in the HH^PT Lagrangian. HH^PT predictions for the ratios in 
Eqs. (0121) are valid for either interpretation, and so cannot distinguish between these two 
scenarios. 

However, if D s0 (2317) (D s i(2460)) is modeled as a bound state of a D (D*) meson 
and a kaon the predictions for the electromagnetic branching ratios will be different. In 
this scenario, instead of adding the even-parity heavy-quark doublet to the Lagrangian by 
hand, the dynamics of the theory containing only the ground state heavy-quark doublet 
and Goldstone bosons generate the observed _D s o(2317) and _D s i(2460). This interpretation 



has been pursued in Refs. |47|, |48|, |49|, |50j, In this scenario the binding energy is only 



about 40 MeV, so the mesons in the hadronic bound state are nonrelativistic. The decay 
rates can be calculated by convolving the unknown nonrelativistic wavefunction with leading 
order HH%PT amplitudes for D^K — »■ D^^y and D^*'K — > D^n . Dependence on the 
bound state wavefunction drops out of the ratios in Eqs. (0|21)- The resulting predictions for 
these branching ratios are much larger than experiment. Furthermore, the branching ratio 
for .D sl (2460) — ► D s / y is predicted to be the smallest of the three, in direct conflict with 
experimental observations. A DK molecular interpretation of the -D s o(2317) and -D s i(2460) 
is disfavored by the existing data on electromagnetic branching fractions. 



II. ELECTROMAGNETIC AND STRONG DECAYS IN HHxPT 

In the heavy-quark limit, hadrons containing a single heavy quark fall into doublets of 
the SU(2) heavy-quark spin symmetry group. Heavy hadrons can be classified by the total 
angular momentum and parity quantum numbers of their light degrees of freedom, j p . The 



ground state doublet has j p = | and therefore the mesons in 
states. In HH^PT, these states are combined into a single field 



the doublet are and 1" 



H a = l±l (fl£ 7 „ - H a75 ) , (4) 
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where a is an SU(3) index. In the charm sector, H a consists of the (D°,D + ,D+) ~ 
(cu, cd, cs) pseudoscalar mesons and H£ are the (D *, D + *, vector mesons. The doublet 
with light degrees of freedom j p = | + consists of mesons whose quantum numbers are + 
and 1 + . These are combined into the field |52j 

S a = ^ - S a ) , (5) 

where the scalar states in the charm sector are S a = D 0a and the axial vectors are S£ = D la . 
The relevant strong interaction terms in the HHxPT chiral Lagrangian are [3,13] 

f 2 f 2 B 
C = — Tr^E^S 1 "] + J — -Tr[m 9 £ + tfm q ] 
8 4 

-Ti\H a w ■ D ba H b ] + Ti[S a (iv ■ D ba - 5 S h 8 a b)S b } 

+ gTr[H a H b fibals] + g'Tr{S a S b fibals] + h(Tr[H a S b fibals] + h.c.) 

-*fTr\H a <r> a 'H tt <T ia ,] + ^-Tt\S a <T^S a a^\ . (6) 



The first line in Eq. (JHJ) is the leading order chiral Lagrangian for the octet of Goldstone 
bosons. / is the octet meson decay constant. The conventions for defining S in terms of 
meson fields, the chiral covariant derivative, D ab , and the axial vector vector field, A^ b , are 
identical to those of Ref. j^. Here m q is the light quark mass matrix. The second line 
of Eq. (Jf)J) contains the kinetic terms for the fields H a and S a and the couplings to two 
and more pions determined by chiral symmetry. The parameter 5,sh is the residual mass 
of the S a field. The H a residual mass can be set to zero by an appropriate definition of 
the H a field, and this convention is adopted here. Then 5,sh is the difference between the 
spin-averaged masses of the even- and odd-parity doublets in the heavy quark limit. The 
third line contains the couplings of H a and S a to the axial vector field A^ b = —d^Ti ab /f + .... 
These terms are responsible for transitions involving a single pion. The couplings g, g', 
and h are parameters that are not determined by the HHxPT symmetries. The last line 
in Eq. (jHJ) contains operators that give rise to l~-0~ and l + -0 + hyperfine splittings, which 
are Ah and A5, respectively. Since the splittings should vanish in the heavy quark limit, 
As ~ Ah ~ A-QCT>/ m Q- The parameters A5 and A# are independent in HH^PT so there is 
no relation between the hyperfine splitting in the even- and odd-parity doublets. In parity 
doubling models, Ah = A$ at tree level, in agreement with the observation that hyperfine 
splittings are equal to within 2 MeV. 
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Electromagnetic effects are incorporated by gauging the U(l) em subgroup of SU(3)l X 
SU(3)r and adding terms to the Lagrangian involving the gauge invariant field strength, 
F^ u . Gauging derivatives in Eq. (jUJ) does not yield terms which can mediate the (0 + , 1 + ) — > 
(0~, 1~) electromagnetic decays at tree level. The leading contribution to these decays comes 
from the operator 

C= e -^TT[H a S b a^]F, u Qt, (7) 



where Q\ a = ±(£Qft + tfQ£) ba , £ 2 = £, and Q = diag(2/3, -1/3, -1/3) is the light quark 
electric charge matrix. A tree level calculation of the decay rates using Eq. ((7J) shows that 

Ik. " 



r[i+- 


- la 7] = 


2 ~ 7 

3 ae q P - 


r[i+- 


- o- 7 ] = 


1 ae 2 J 2 - 
3 q 1 






9 ~ m 


r[o+- 


- Ki\ = 


ae 2 (3 2 
H m, 



n 1 



la II, 13 



3 



k 7 r, (8) 

y 'a 

where e q is the electric charge of the light valence quark, a is the fine-structure constant, 
and j3 is the unknown parameter in Eq. (J7J). The three-momentum of the photon in the 
decay is k 7 and mjp is the mass of the heavy meson with quantum numbers J^. In the 
heavy-quark limit, the members of each doublet are degenerate and the phase space is the 
same for all three decays. If differences in phase space are neglected the decay rate ratios 
are T[l+ — > 1^7] : T[l+ — > 0~j] : r[0+ — > 1~7] =2:1:3. Differences in the phase space 
factors are formally 0(l/mg) but in practice it is critical to include these effects to make 
sensible predictions. For the charmed strange mesons using the physical masses gives 

r[Aa(2460) D* s j] = GeV) 2 15.6 keV 
r[ J D sl (2460) -> D sl ] = GeV) 2 18.7keV 

r[D a0 (2317) -»• -D*7] = (/3GeV) 2 5.6keV. (9) 

The rates are then in the following ratios: 

r[£> sl (2460) -> D* sl ] : T[D sl (2A60) -> D s -f] : r[D s0 (2317) -> D*j] = 0.83 : 1.0 : 0.30 . (10) 

Note that the rate for Z) sl (2460) — ► D s ^y, smallest in the exact heavy-quark limit, is actually 
the largest when phase space effects are included since |k 7 | is largest for this decay. 
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To compare with the ratio measured by Belle, the isospin violating strong decays must 
be calculated. These decays proceed through 77 — 7r° mixing. The result is [14] 



2,-Kp m Dsl(2 460) 

17.0 keV if f = f w = 130 MeV 
9.8 keV if / = /„ = 171MeV 



h 



r[D s0 (2317) - D S A = ff 2 mDa Egolp^l 

9 f 16.9 keV if f = U= 130 MeV 
= fc 2 4 , (11) 

[ 9.8 keV if / = /„ = 171MeV 

where 6 1 = y/3/2(md — m u )/(2m s — — m u ) = 0.01 is the r\ — n° mixing angle. E n o and 
p n o are the energy and three-momentum of the 7r°, respectively. At tree level / = f n — f v . 
The difference between the two predictions provides an estimate of the uncertainty due to 
higher order £7/(3) violating effects. 

The branching fraction ratios in Eqs. $I]|2J) depend only on the ratio (3 2 /h 2 at leading 
order in HHxPT. To obtain h 2 separately, a measurement of an excited strong decay width 
is needed. Currently h 2 cannot be extracted from the strange sector because only loose 
experimental bounds on r[D s o(2317)] and r[D s i(2460)] exist. Until measurements of these 
widths are dramatically improved, h 2 can be estimated using data on nonstrange even-parity 



D meson widths. CLEO 55| has observed preliminary evidence for the D® (J p = 1 + ) meson. 
(Here the superscript refers to the particle charge.) More recently, Belle has reported 
observing; even- parity D° ( J p = 0+) and D® states. Finally, FOCUS has observed broad 
structures in excess of background in the D + 7r~ and D°tc + invariant mass spectra. FOCUS 
does not claim to observe an excited charm resonance but does fit the excess with a Breit- 
Wigner to determine the resonance properties required to explain their data. The masses 
and widths reported by all three experiments are collected in Table |l] The experiments all 
quote several errors which have been combined in quadrature for simplicity. Note that the 
CLEO and Belle measurements of the are consistent with each other while the central 
value of the mass obtained by FOCUS is 99 MeV higher than the central value of the 
Belle measurement. Furthermore, the FOCUS D® mass is actually greater than the mass of 
the .D s o(2317). If the effect observed by FOCUS is a scalar D resonance, it seems implausible 
that this resonance is related to the D s0 (2317) by £7/(3) symmetry. Therefore, the FOCUS 
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Experiment 


Particle(J p ) 


Mass (MeV) 


Width (MeV) 


CLEO [55] 




246ll| 


290l«£° 


Belle [56] 


D§(0+) 


2308 ± 36 


276 ± 66 






2427 ± 36 


384±gg 


FOCUS [57] 


D °(0+) 


2407 ± 41 


240 ± 81 




^o + (0 + ) 


2403 ± 38 


283 ± 42 



TABLE I: Masses and widths of even-parity non-strange charmed mesons, -Dj , where Q is the 
electric charge. 

data will not be used to estimate h 2 . Even the masses obtained by CLEO and Belle are 
large compared to expectations based on 577(3) symmetry. Combining the strange sector 
+ -0~ and 1 + -1~ mass splittings with SU(3) symmetry leads to the prediction that the Dq 
mass is 2212 MeV and the D® mass is 2355 MeV |l3|. 

Applying the leading order expression for the decay widths of the nonstrange + and 1 + 
mesons 



2nf 2 

h 2 



m D +* 



1 m£,ot 



2 m D o 



m D+ E 2 



2ixf 2 \m D o 



\Ptt- 



1 m D o 2 



2m n o 



(12) 



to the CLEO and Belle data yields h 2 = 0.39 ±0.13 from the 0+ decays and h 2 = 0.49 ±0.14 
from the 1 + decay. The error in each case is obtained by adding in quadrature the uncertainty 
in the decay rate from varying the mass within the allowed range and the experimental error 
in the decay rate. If the two results are averaged h 2 = 0.44 ± 0.11. This estimate of h 2 is 
consistent with the bound h 2 < 0.86 extracted from an Adler-Weisberger type sum rule for 



7rB scattering 



581 ] . (To obtain this bound g = 0.27 |5J] is used in the result of Ref. [58].) It 
is also consistent with a calculation of h = —0.52 ±0.17 obtained using QCD sum rules in 
Ref. Q. 

The lowest order HH%PT prediction for the ratio measured by the Belle collaboration is 



Br(^ 1 (2460)^D s7 ) ( PGeT 
Br(L> s i(2460) -> D* s 7i°) \ h 



x 



1.1 
1.9 



if / = /* 
if / = /„ 



130 MeV 
171 MeV 



(13) 
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Averaging the results of the two Belle measurements, (j3 GeV/h) 2 = 0.40 ±0.08 (0.23 ±0.05) 
or \(3\ = 0.42 ± 0.07 (0.32 ± 0.05) GeV -1 , where / = /„. (/ = f v ). The error is estimated by 
first combining the statistical and systematic errors in quadrature for each measurement, 
then combining the two measurements assuming they are independent. The extracted values 
for j3 and h are consistent with expectations based on naturalness. Plugging the value of 
(f3GeV/h) 2 into expressions for the unobserved electromagnetic decays yields the following 
predictions for the branching fraction ratios: 

6^(2460)^) Br( C , (2317) 

Br(L> sl (2460) -> D*n°) Br(D fl0 (2317) -> D s n°) v ; 

Both predictions are in excess of bounds quoted by the CLEO experiment. Heavy-quark 
spin symmetry predicts branching ratios for the electromagnetic decays _D s i(2460) — > _D*7 
and D s o(2317) — > D*^ that are more than a factor of two in excess of the experimental 
upper limits. 

In the rest of this section the leading corrections to both electromagnetic and strong 
decays are analyzed. Because Aqcd/^c ~ 1/3 corrections to heavy-quark spin symmetry 
predictions can be rather large for charm hadrons. These corrections can be systemati- 
cally analyzed using HH%PT. For example, the pattern of deviations from heavy-quark spin 
symmetry predictions for the one-pion decays of excited D-wave charm mesons |5j| can be 
understood by analyzing the structure of spin-symmetry violating operators appearing at 
0(l/m c ) in the HHxPT Lagrangian 6^|. A ratio of decay widths for which the 0(l/m c ) 
correction vanishes agrees well with data. There is another ratio for which the leading order 
heavy-quark spin symmetry prediction fails rather badly. In this case the leading 0(l/m c ) 
correction is multiplied by a large numerical coefficient. Thus HHxPT is a useful tool for 
determining the robustness of predictions based on heavy-quark spin symmetry. 

Spin-symmetry violating operators that contribute to 5* — > H transitions must have 
the Dirac structure TrfHV^S^] or Tr \Ho^ u S^ a ] . Operators with HS conserve spin sym- 
metry, while operators with Ha^ u S violate spin symmetry. Operators with H'y^S and 
#y 75 S are redundant since H^S = H\{^,i)}S = v^HS and Ha^S = ~H\{a^^}S = 
—e^ ua/ 3V a H'y f3 'y 5 S, while H^ 5 S = 0. Spin-symmetry violating operators will be of the form 
Tr^HcrV ST]. T must be 7° or 75 since the trace vanishes for F — 1, while T = 7^75 and a af3 
are redundant because 

Tr[i7a^S 7 Q 7 5 ] = v a Tr[Ha^ S-f 5 ] 
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Reparametrization invariance 61] forbids operators with derivatives acting on the H or S 
fields j^. For S — > decays, the lowest dimension, parity conserving, spin- symmetry 
violating operators are 

C = % -^-Tr\H a <T^Sa7s]F^e^ a p + ^^Tr[tf a a^ a7 Q ^ a i^ + h.c. . (15) 
8m Q 8m Q 

The 1 / itlq dependence (expected for any operator which violates heavy-quark spin symme- 
try) is explicit. The factors of i are required by time reversal invariance. The first operator 
in Eq. (|15j) and the leading operator in Eq. (JJJ) have mass dimension 5. The second operator 
in Eq. (|15|) has mass dimension 6, so (3" has mass dimension —1 and is expected to scale 
like 1/A X ~ GeV -1 . Since 2a^ u d^F ai/ = a^daF^ for an abelian field strength there is a 
unique way of contracting indices in this operator. Note that there is a unique dimension 6, 
spin-symmetry conserving operator proportional to Tr^HaSba^Q^iv ■ dF^ u . This operator 
gives slight deviations from the ratios in Eq. fllOl) since its contribution is suppressed by 
|k 7 |/A x and |k 7 | differs for the three decays due to hyperfine splittings. These corrections 
should be smaller than corrections coming from operators in Eq. (|15|) so they are neglected 
in what follows. 

Power counting is used to determine the importance of higher dimension operators in 
the Lagrangian. HHxPT is a double expansion in Aqcd/^-q an d Q/A x , where Q ~ p ~ 
ma ~ rriK- Two additional mass scales appearing in the Lagrangian of Eq. (jHJ) are the mass 
splitting between the H and S doublet fields, 5sh, and the hyperfine splittings within each 
doublet. In the heavy-quark limit, the S field propagator is proportional to 

2(vk-6 SH )- (16) 
If Ssh were to scale as Q°, then the S propagator could be expanded in powers of v ■ k/SsH 
since v ■ k ~ Q. In the strange sector loops receive important contributions from momenta 
~ rrix = 495 MeV. Numerically, 5sh ~ 350 MeV in the strange quark sector and ~ 430 MeV 
in the nonstrange sector, so expanding in v ■ k/5sn is a poor approximation. Therefore, 
$sh ~ Q is required. The hyperfine splittings are also treated as ~ Q since numerically 
these splittings are ~ 140 MeV which is ~ m n . 

There are SU(3) violating corrections to the decay rates from operators such as 
Tr[H aSba^F^Qljn^. These operators will give the same correction to all three elec- 
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tromagnetic decays in Eqs. (0121), so their effect can be absorbed into the definition of /3. 
However, if one were interested in relating the electromagnetic decays of strange and non- 
strange heavy mesons these operators must be included explicitly. 

The leading operator in Eq. (J7J) is order Q because of the derivative in F^ u . The first 
operator in Eq. (fTK|) is treated as ~ Q 2 because it is suppressed by Aqcd/^q relative to 
the leading operator. The second operator in Eq. (fTKJ) has two derivatives and is also I/tuq 
suppressed. It is treated as ~ Q 2 . The correctness of this power counting is confirmed by 
the calculation of one-loop corrections to the decay, since in order to properly renormalize 
these diagrams both counterterms in Eq. (fTo'j) are needed. In loop diagrams, integrals scale 
as Q 4 , the propagators of the H and S fields as ~ Q^ 1 and the propagators of Goldstone 
bosons as ~ Q~ 2 . The leading couplings of the H and S fields to kaons and pions are ~ Q 
and the couplings of the photon to the kaons and pions are ~ Q. Calculations of the loop 
corrections are performed in v ■ A = gauge where the leading coupling of the photon to 
the heavy meson fields vanishes. Finally, there is an ~ Q° coupling of two heavy meson 
fields to a Goldstone boson and photon which comes from gauging the derivative couplings 
of the heavy meson fields to Goldstone bosons. The HSttj vertex obtained by gauging the 
derivative on the pion field in the leading HSn coupling vanishes in v ■ A = gauge. With 
these power counting rules, the loop diagrams shown in Fig.l give an 0(Q 2 ) contribution to 
the S — > if 7 decays. Double lines are S fields, solid lines are H fields and the dotted lines 
are Goldstone bosons. For D s decays the Goldstone bosons in these loops are K + and the 
virtual heavy mesons are neutral .D's. 

Including both the loop diagrams and tree level insertions of the operators in Eq. |T5|) 
yields: 

f[l a + ^l-7] = l + ^(-^+F\m lt -m 1 -,m lt -m lt AKlM,fi\] (17) 



f[l+^0- 7 ] = 1 + A ( e ^L + e QP"\ k f\ + F \ m _ m m _ m |kJ M, ^ 

a e q (3 \m Q 2m Q L 1 » h % 7 " 

r[0+ — > = l + -^(^^ — GQ ^ ^ 7 + F |m n + — m n - ,m n + - m,+ , IkJ, M, u 
e q f3\m Q 2m Q L °« o b o a i b i ti> ,h>^ 

Here f = y nlo /T LO , where the T LO are given in Eq. (JHJ). The SU(3) index a refers to 
the external heavy mesons while the index b refers to the mesons inside the loop. M is 
the mass of the virtual Goldstone boson. For heavy-strange decays the external particles 
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FIG. 1: One-loop chiral corrections to the electromagnetic decays S — > Hj in v • A = gauge. 
Double lines are S mesons, solid lines are H mesons, dashed lines are Goldstone bosons and the 
wavy line is the photon. 

are heavy-strange mesons; a = 3, the Goldstone boson is a K + , and the heavy- mesons 
inside the loops are neutral heavy mesons with 6=1. f is expanded to 0(Q). The 
function F[Ai, A 2 , |k 7 |,M, fj] is given in the Appendix. The loop graphs are regulated in 
dimensional regularization, counterterms are defined in the MS scheme and the dimensional 
regularization parameter is \i. All \i dependence is canceled by the implicit \x dependence of 
the renormalized couplings /?, /?', and ft". 

An NLO calculation of the electromagnetic branching ratios also requires 0(l/m c ) cor- 
rections to the decays T[D s0 {2317) -> D s n°] and r[L> sl (2460) -> D*ir ]. The leading spin- 
symmetry violating operator contributing to these decays is 

£ = ^Tr[H a a^S b r\Al . (18) 

Because of the 1/toq suppression this operator is considered 0(Q 2 ). The one-loop diagrams 
contributing to S — > Hn transitions are subleading at 0(Q 3 ). The decay rates to NLO are 

r[0 3 + - 3 -7r°] = (h + 3—) ^-^E^M. (19) 
V m Qj 37r/ 2 m + 

Earlier in this section the data from Dq and D® decays was averaged to extract h 2 . Including 
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the leading correction it is possible to fit h and hi separately, extracting h = 0.69 ±0.09 and 
h'/m c = -0.019 ±0.034. 

The NLO expression for the branching fraction ratios of heavy-strange mesons is obtained 
by combining Eq. (J 19)) with Eq. (|17|) and Eq. (jHJ). The result is 



Br[l 3 + ^l 3 - 7 ] 2tt«/ 2 /3 2 Ik, 



BrfK 



IsTrO] 



9 2 /i 2 E^o 



x 



BrfK 



2h' 6e Q 0' 6 _ 
l + - 7 + ^L-^F 

hm Q (3m Q 

-> 3 - 7 ] naf 2 2 m 



m,+ — to,- , m,+ — m,+ , IkJ, m^+, u 



Br[K 



IsTrO 



9W m x -El oPn o 
2h! 6e Q 0' 3e Q ^"|k 7 | 6 



hm Q (3m Q 0m Q 



m x + — m x - , 771^+ — m + , |k 7 |, m K +, li 



Br[0+ -> I3 7] , Traf^m,- |k 7 | 



Br[0 3 : 



OgTT ] 



30 2 /i 2 m 07 E 2 oPw o 



(20) 



6/i' 6e Q /3' 3e Q /3"|k 7 | 6 



+ 



- —F 





m n + — m n - , m n + — m,i , IkJ, TOir+, W 



V (3m Q 0m Q 

where e q has been set to e s = —1/3. Applying the formulae in Eq. (|20|) to the experimentally 
observed ratios gives 

Br(D sl (2460)^D; 7 ) _ 5g 0* 
Br(L> sl (2460) -> D s *vr°) ' /i 2 

/ 1.43 /i' 2.86/3' 0.18 g/i 2.94 + gIg ^ 7 /i\ 
1 + — : — + — J- + =F ~ - I < 0.16 



h 















Br(D sl (2460) -> D sl ) 



2 

1.90 



Br(D s i(2460) -> 

1.43 /i' 2.86/3' 0.63/3" 0.03 ^/i 



1 + 



2.40+0.73 „/ 



-0.59 



/3 















0.44 ±0.09 



Br(^o(2317) -> £ s7 ) n . 7 ^! x 
Br(D s0 (2317) -> D S 7T°) ' h 2 



4.29 h' 2.86 0' 0.28/3" 0.37 gh 3.81±% f 5 g'h 



1 



/1 



+ 







< 0.059. (21) 





Here /3 and /3" are measured in units of GeV -1 and hi in units of GeV. All other quantities are 
dimensionless. The charm quark mass is m c = 1.4 GeV, the renormalization scale is li = 1 
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GeV, and cq = e c = 2/3. For the loop corrections with kaons the meson decay constant is 
f = fx, while for the strong decays / = f v . (If f n is used in the strong decays, then the 
branching fraction ratios in Eq. (|2U|) should be multiplied by f%/f% = 0.58.) The masses 
used for the virtual nonstrange even-parity heavy mesons in the loops are m + = 2308 ± 36 
MeV and mj+ = 2438 ±29 MeV, where the first number is the nonstrange + mass measured 
by Belle and the second is the average of the nonstrange 1 + mass measured by CLEO and 
Belle. The uncertainty in the coefficient of g'h//3 in Eq. (j21j) is due to the uncertainty in 
the masses of the -Dq and D®. 

The result depends on seven parameters: g, g' , h, h' , j3, j3' , and (3". The coupling q i s 
constrained to be 0.271q3 from a next-to-leading order HHxPT analysis of D* decays [54]. 
h and h! are extracted from the nonstrange decays, leaving four unknown parameters. Since 
there are only three constraints coming from experiment, further analysis requires additional 
assumptions to constrain the parameter space. 

To illustrate how the current data is consistent with natural size parameters the follow- 
ing situation is considered. The contribution from (3" is neglected since in Eq. (|2*T|) (3" is 
multiplied by a coefficient that is much smaller than the coefficients multiplying h' and (3' . 
(The smallness of this coefficient is due to the factor |k 7 |/m c .) g, h, h', and the branching 
fraction ratio measured by Belle are set to their central values: 0.27, 0.69, -0.019 m c , and 
0.44, respectively. Ranges for the remaining parameters (/?, (3', and g') are extracted by 
varying the branching ratios in Eq. (J2J) between and their upper limits. There are two 
solutions since the formulae for the electromagnetic decay rate is quadratic in (3. The results 
are 

0.70 < GeV < 0.86 - 0.01 < (3' < 0.01 0.32 < g' < 0.40 
-0.62 < (3 GeV < -0.46 - 0.01 < (3' < 0.02 - 0.25 < g' < -0.16 . (22) 

Note that the ranges quoted in Eq. (|2*2*|) do not include errors due to the uncertainties 
in the parameters g, h, and h' or the masses of the Dq and D®. h! is highly uncertain 
because of the uncertainty in the masses and widths of the Dq and D® used to extract it. 
The loop contribution proportional to g'h/(3 is also sensitive to the masses of the Dq and 
D\ that appear as intermediate states. The ranges given in Eq. (|2*2*j) do not reflect these 
uncertainties and do not exhaust the possible parameter space. Instead, they are simply 
illustrative of natural size parameters consistent with existing data. 
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When more data on excited heavy meson systems becomes available, the formulae in 
Eq. (j2*Uj) could be used to test models that make predictions for the parameters in HHxPT. 
In parity doubling models, g = —g' and h — 1 at tree level Q|. The authors of Ref. 
note that h can be renormalized away from its tree level value and allow this parameter 
to vary in their analysis of strong decays. The tree level result h = 1 exceeds the value 
extracted from excited nonstrange decays in HHxPT. Another theoretical framework which 
makes similar predictions for the coupling constants g, g', and h is the algebraic realization of 
chiral symmetry j^jj]. Applying this theory to heavy mesons leads to the predictions 

g' = —g and g 2 + h 2 = 1. Using g = 0.27 in this relation gives h 2 = 0.93 which is also 
larger than extracted from Eq. (|19|) . While the predictions for h are not in agreement with 
available data, the condition g = —g' is consistent with available data but not required. 

Eventually the even parity B s states will be observed and all electromagnetic branching 
fractions for heavy-strange mesons will be measured. Then the parameter space will be over- 
constrained and HHxPT for excited heavy mesons can be tested decisively. Furthermore, 
the extracted values for g, g', and h can be compared with predictions from parity doubling 
models and algebraic realizations of chiral symmetry. At the present time, observed viola- 
tions of leading heavy-quark spin symmetry predictions are consistent with what is expected 
from loop effects and higher order operators appearing in the HHxPT Lagrangian. 



III. ELECTROMAGNETIC DECAYS AND D K MOLECULES 

The unexpectedly low masses of the D s0 (2317) and _D s i(2460) have prompted speculation 
that these states are unconventional. Two common proposals are that these mesons are 
csqq tetraquarks or hadronic bound states of D and K mesons. This section addresses 
the question of what the decays reveal about the internal structure of the _D s o(2317) and 
_D s i(2460). In the analysis of the previous section the only information about the states 
needed to construct the HH;\PT Lagrangian is the assumed SU(3) and j p quantum numbers 
of the light degrees of freedom in the hadrons. A constituent quark in a P-wave or an exotic 
with two light quarks and an antiquark both have j p = \ ■ Both states are represented by 
a field like that in Eq. (JHJ). Analysis of electromagnetic and strong decays within HIPyPT 
is identical for both states, though the coupling constants f3, /3', h, etc. would be different 
for the two states. Since these coupling constants are unknown in either case, the HHxPT 
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predictions for electromagnetic and strong decays cannot distinguish between exotic csqq 
and conventional cs P-wave states. Of course, if the Z) s o(2317) and ^^1(2460) are csqq 
states then in the quark model there should be distinct cs P-wave mesons with the same 
quantum numbers. These states could be very hard to detect, however, if they are above the 



28, 



3- 



DK threshold. Mixing between the conventional and exotic mesons is also likely 

However, if the D s o(2317) (D s i(2460)) is a bound state of and K mesons then the 
HHxPT predictions for electromagnetic and strong decays will be different. For a hadronic 
bound state of a D or D* and a kaon, one could in principle calculate the bound state masses 
and other properties from the HH^PT Lagrangian with the field H a alone. There have been 
attempts to g enerate the D s0 (2317) and D sl (2A60) as resonances in a unitarized meson 
model as well as by solving Bethe-Salpeter equations in relativistic, unitarized chiral 



perturbation theory 



50 



5l|. Producing a bound state requires resumming an infinite 



number of Feynman graphs in HH;\RT and the renormalization of these graphs requires 
introducing higher order operators whose renormalized coefficients are unknown. Such a 
calculation will not be attempted in this paper. Instead the DK molecular picture will 
be tested by simply assuming that strong forces between and K mesons give rise to 
the _D s o(2317) and £) s i(2460) and determining what this implies for the decay rates. If the 
_D s o(2317) (£) s i(2460)) are bound states of D^*'K then the characteristic momentum of the 



constituents is p ~ y/2pJ3 « 190 MeV, where \x is the reduced mass and B is the binding 
energy. The DK molecule can then be modeled as a nonrelativistic bound state since 
relativistic corrections are suppressed by v 2 = p 2 /M] i = 0.15. Strong and electromagnetic 
decays can be calculated in terms of the unknown bound state wavefunction. Even without 
any knowledge of these wavefunctions it is possible to make predictions for the decay ratios 
in Eqs. It turns out that these predictions disagree with data so interpreting D s o(2317) 

and _D s i(2460) as DK molecules is disfavored. 

Ref. 23] advocates using the radiative decays of the _D s o(2317) and -D s i(2460) to deter- 
mine the nature of these states and calculates radiative and strong decays within a non- 
relativistic quark model. Predictions for the branching fraction ratios in Eqs. ((TJEJ) are in 
the same proportion as leading order heavy-quark symmetry predictions, though they are 
approximately 45% larger than the leading order predictions obtained in section II. The 
quark model expectation for the total widths of the _D s o(2317) and D s i(2460) is O(10 keV), 
consistent with Eq. pijl . However, the conclusion of Ref. [2J| states that a D^K molecule 



17 




FIG. 2: Leading order diagrams for D^K bound states decaying into Dg 7. The shaded oval 
represents the D^K bound state wavefunction. 

should have a width of 0(1 MeV) and that the electromagnetic transitions should be absent. 
The analysis that follows is consistent with the first conclusion but not the second. Below 
it is demonstrated that the electromagnetic branching ratios of a D^*'K molecule are large 
and are in worse agreement with experiment than the nonrelativistic quark model. 

The _D s o(2317) (D s i(2460)) is assumed to be an S"-wave 1 = bound state of and K 
mesons. The matrix elements for the electromagnetic decays of the _D s o(2317) and _D sl (2460) 
are given by: 

M[D sQ {2m) -> D* sl \ = .P?— f ^-J DK ( p )M[D(p)K(-p) -> D sl \ 

V ™d s0 J (27T) 3 

M[D sl (24Q0) ^ Di*h] = ^/^^(pl^WN-^i (23) 

V rn Dsl J (2vrf 

Here p is the three-momentum of the D^*' meson in the bound state and ^wjf(p) is 
the bound state momentum-space wavefunction. Although calculation of the bound state 
wavefunction is nonperturbative, the typical momentum is small enough that the amplitudes 
M. [D^K — > D^'j] are perturbatively calculable in HHxPT. The leading order diagrams for 
the decay rates in Eq. (J23j) are shown in Fig.|21 The shaded oval on the left hand side of these 
Feynman diagrams represents the D^*'K molecule. In Figs. |2^-b the dashed line is a K + , 
and the vertex involving the photon comes from gauging the D*-D s -K + coupling (Fig. |2K) 
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or the K + kinetic term (Fig. |2p). in Figs.|2fc-d, the photon coupling comes from gauging the 
heavy meson kinetic term. There are also diagrams like Figs. Et-d where the photon heavy- 
meson coupling comes from a term in the Lagrangian proportional to TT[HbH a cr^ u Q^ ab ] F^ u , 
but these only contribute in the P-wave channel. 

The graphs in Fig.|2t-d are nonvanishing in the K°D + channel, but are equal and opposite 
in sign so the contribution in this channel vanishes. The graph in Fig. |2fc vanishes in the 
K + D° channel. The amplitudes are 

2eg 



y/m D *m D * 



M[D*{p)K{-p) - D* sl ] 
M[D*(p)K(-p) - D 8l ] 
M[D{p)K{-p) D* sl ] = jm^> 



Pk(PK-PjY +g »5_ PK vfl 



y/m D *m D .— 

2eg_ 
f 



f L PK-P-i 

PxiPK-p^y 

PK ■ P-y 

Pk(Pk-Pj) 1 



v-pk 



v^p u K 

V-Pk 

v»p v K 



(24) 



PK-P-/ V-p K _ 

Here px and p 7 are the kaon and photon four-momentum, respectively. The polarization 
vectors for the photon, D*, and D* are denoted e, ei, and e 3 , respectively. It is easy to 
check that the amplitudes respect the QED Ward identity. These expressions are inserted 
into Eq. (j2HJ), p^ is set to Ekv^ + p M , p M = (0, — p), and the matrix element is expanded to 
lowest order in p. Because of the rotational symmetry of the S'-wave wavefunction, ip(p), 
terms linear in p vanish. There are corrections to the amplitudes from higher orders in chiral 
perturbation theory that are 0(m 2 K /A^ c ) and relativistic corrections of 0(v 2 ). The errors in 
the predictions for the decay rates could be as large as 50%. 
The results for the decay rates are 

r[D sl (2460) - W ^(m^m Dt \ 
r[Aa(2460) -+ D sl ] 
r[ J D s0 (2317) -> D* sl ] 



3/ 2 


V m ki 


Ag 2 a 


( m D o*m Ds 


3/ 2 


V m ki 


Ag 2 a 


( m D om D *^ 


f 2 


V m ko / 



) \lpD*K\ 

I |^(0)| 2 |k 7 | 
|^ DX (0)| 2 |k 7 |. 



(25) 



Here i/)dk(0)[ 1 Pd*k(0)] is the wavefunction at the origin for the D sQ (2317) [D 8 \ (2460)]. In the 
heavy-quark limit the partial width ratios are again r[.D s i(2460) — > D*^] : r[.D s i(2460) — > 
Dg^] : r[D s0 (2317) — ► D*j] = 2:1:3. However, the decay rates are proportional to |k 7 | 
instead of |k 7 | 3 . This important difference in the kinematic factors leads to a very different 
prediction for the relative sizes of the partial widths than obtained in Eq. ()10j). In this case 

r[L> sl (2460) -> D*j] : r[L> sl (2460) -> D s j] : T[D s0 (2317) -> D* sl ) = 1.57 : 1 : i^l.58, 
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FIG. 3: Leading order diagram for D^K bound states decaying into D^ir . The dashed line 
from the bound state is a K, the dashed line in the final state is an rj which mixes into a 7r°. 

where = \iPdk(Q)\ 2 /\iPd*k(0)\ 2 is expected to be « 1. In this scenario T[D sl (2A60) — > 
Dg^] is the smallest decay rate rather than the largest. 

To compare with the measured branching ratios the strong decays must also be calculated. 
The leading order diagrams are shown in Fig. El All three diagrams depict D^K — > D^rj 
followed by rj — 7r° mixing, which is represented by a cross on the dashed line in the final 
state. The vertex for the graph in Fig.EK comes from the chirally covariant derivative in the 
heavy meson kinetic term. This graph contributes in the S-wave channel while Figs. H2>c 
contribute to the P-wave channel only. The results for the decay rates are 



r[ J D sl (2460) -> D s ti ] = y ^ ) \ip D * K {0)\ |iv| 



rm / oq17 n n 01 3(m K + E^ye 2 (m D m Ds \ 2 
r D a0 (2317) -> D s ii>\ = — 3 \iPdk{0)\ Ptt" (26) 



If / = f n = 130 MeV then the bounds r[D sl (2460)] < 7 MeV and r[Z> a0 (2317)] < 7 
MeV imply \ip DMK (0)\ 2 < (52MeV) 3 . If instead / = /„ = 171 MeV then \ip DMK (0)\ 2 < 
(75 MeV) 3 . In either case the bounds on the wavefunctions are somewhat smaller than 
expected: IV'dMa'WI 2 ~ |p| 3 ~ (190 MeV) 3 . Since this is only an order of magnitude 
estimate, the bounds on \ip d(*) k{^)\ 2 are n °t a problem for the DK molecular interpretation. 
However, they do imply that if the _D s o(2317) and _D s i(2460) are D^K molecules the states 



23, 



27]. 



should not be much narrower than the present upper limits 

Since the wavefunction squared cancels in the ratio of strong and electromagnetic decays 
the electromagnetic branching fractions can be predicted: 

Br(D sl (2460) -> D* sl ) 



Br(P sl (2460) -> D*ir°) 



3.23 
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Br(/J sl (2460) - D sl ) _ ^ 



Br(D s i(2460) -> D*tt ] 
Br(/J s0 (2317) -> £.7) 



2.96. (27) 



Br(A, (2317) -> D s ir°) 

In this calculation a = 1/137, 9 = 0.01, g = 0.27, = 159 MeV in the electromagnetic 
decays and f = f v = 171 MeV in the strong decays. If instead / = f n = 130 MeV is used 
in the strong decays the predicted branching fraction ratios are smaller by a factor of three. 
While the branching fraction ratios are quite sensitive to the choice of /, in any case they are 
much too large compared to experiment. Also, the relative sizes of the branching fraction 
ratios are in disagreement with experiment, since the second branching fraction ratio in 
Eq. (j2Sj) is predicted to be smallest, not largest. Note that the possibility of these states 
being mixtures of quark level bound states and DK molecules 3, 0] is also disfavored since 
this would enhance the first and third ratios in Eq. (|2T)j) relative to the second, whereas in 
reality these ratios are suppressed relative to the leading order prediction in Eq. (fTUj) . 



IV. CONCLUSIONS 



In this paper, corrections to electromagnetic and strong decays of -D s o(2317) and 
-D.,1 (2460) are calculated in HHxPT. The corrections depend on a number of unknown 
or poorly determined coupling constants. These predictions can be consistent with Belle 
and CLEO data with coupling constants of natural size. Serious tests of the HH%PT de- 
scription of the _D s0 (2317) and -D s i(2460) will require more data on the electromagnetic 
branching ratios of the even-parity charmed strange mesons and the strong decays of their 
nonstrange partners as well as the decays of even-parity bottom strange mesons yet to be 
observed. The work in this paper provides further stimulus for better experimental measure- 
ments of charmed strange decays as well as discovery of their bottom strange counterparts. 
Once better data becomes available, it would be interesting to test models of chiral sym- 
metry breaking which make specific predictions for the coupling constants appearing in the 
HH^PT lagrangian. 

This paper also tests the hypothesis that the _D s o(2317) and -D s i(2460) are molecular 
bound states of DK and D*K molecules, respectively. In this scenario, these states are 
sufficiently nonrelativistic that HH^PT can be used to predict the decay rates at lowest 
order. Furthermore, bound state wavefunctions cancel out of predictions for the observed 
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branching fraction ratios so absolute predictions can be made. These predictions are in much 
worse agreement with data than leading order HH;\PT predictions. Specifically, predictions 
for all the branching fraction ratios are larger than observed and the branching fraction 
for the only observed electromagnetic decay is predicted to be the smallest of the three 
possible decays rather than the largest. Therefore, a molecular interpretation of these states 
is disfavored by available data on electromagnetic decays. 

R.P.S. and T.M. are supported in part by DOE grant DE-FG02-96ER40945. T.M. is 
also supported in part by DOE grant DE-AC05-84ER40150. T.M. would like to thank the 
Aspen Center for Physics where part of this work was completed. 



V. APPENDIX 



The function F[Ai, A 2 , |k 7 |,M, //] is 

gh 



F [Ai, A 2 , |k 7 |,M, n] 



8tt 2 / 2 
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(28) 
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The functions Fi j2 (x) are given by: 
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